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In the framework of open quantum systems, we study the dynamics of a static polarizable two-level 
atom interacting with a bath of fluctuating vacuum electromagnetic field and explore under which 
conditions the coherence of the open quantum system is unaffected by the environment totally. For 
both a single-qubit and two-qubit systems, we find that the quantum coherence can not be protected 
from noise when the atom interacts with a non-boundary electromagnetic field. However, with the 
presence of a boundary, the dynamical conditions for the insusceptible of quantum coherence are 
fulfilled only when the atom is close to the boundary and is transversely polarizable. Otherwise, the 
quantum coherence can only be protected in some degree in other polarizable direction. 
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I. INTRODUCTION 

Quantum coherence is an important concept in quan¬ 
tum theory. The superposition principle of quantum 
states summarize the characteristic features of quan¬ 
tum systems from their classical counterparts, which 
marks the coherence is a fundamental aspect of quan¬ 
tum physics [I| . As the technological progress in the last 
few decades, quantum coherence becomes a powerful re¬ 
source in various aspects, such as low-temperature ther¬ 
modynamics |2 hct, quantum metrology [t], l8| and solid- 
state physics [Sl, [lOj. Besides, the interference phenom¬ 
ena plays a key role in quantum optics experiments [n^ 
in, as well as a series of quantum information pro¬ 
cessing tasks [l5l| . Conventionally, quantum coherence 
is regarded as a physical resource which is associated 
with the capability of a quantum state to exhibit quan¬ 
tum interference phenomena. Recently, Baumgratz et 
al. proposed some rigorous measurements to quantify 
coherence, such as li norm and relative entropy of co¬ 
herence |16l| , which are distance-based measurements for 
such resources. 

On the other hand, a realistic quantum system should 
be regarded as an open system due to the interaction be¬ 
tween the system and its surrounding environment. For 
example, a two-level atom bathing in a vacuum fluctu¬ 
ations of electromagnetic field would suffers from deco¬ 
herence, in which quantum coherence would reduced due 
to the interaction between the atom and the environ¬ 
ment. Considering that the coherence is not only a fun¬ 
damental aspect of quantum physics but also a crucial re¬ 
source in quantum information technology, it is essential 
to find strategies to protect it. Most recently, Bromley et 
al. found that under some suitable dynamic conditions. 
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the quantum coherence is efficiently preserved from noise 
during the entire evolution El- 


In this paper we consider a static two-level atom which 
couples with the bath of fluctuating vacuum electromag¬ 
netic field. Our goal is to protect the quantum coherence 
by setting a reflecting boundary in the electromagnetic 
field, which would changes the vacuum fluctuations. We 
first calculate the evolution of the li norm and relative 
entropy of coherence for a single-qubit system. Two dif¬ 
ferent cases, the electromagnetic field without or with the 
presence of a reflecting boundary, are considered. We are 
going to analyze under which conditions the coherence is 
totally unaffected by the vacuum fluctuations. We find 
that in an unbounded space, the quantum coherence can 
not be frozen during the whole evolution, which is due 
to the interaction between the two-level atom and the 
electromagnetic field. However, with the presence of a 
boundary, the quantum coherence can be frozen when 
the atom is close to the boundary and is transversely 
polarizable, which means that the quantum coherence 
can be protected efficiently. We then study the quantum 
coherence of a two-qubit system with maximally mixed 
marginals [isl . IT^ and find that the result is similar to 
that of the single-qubit case. 


The outline of the paper is as follows. In section II, 
we briefly introduce the methods to quantify coherence 
and discuss the dynamics of open quantum systems for 
a static polarizable two-level atom coupling a bath of 
fluctuating vacuum electromagnetic field. In section III, 
we calculate the evolution of the li norm and relative 
entropy of coherence for a single-qubit system. In section 
IV, we study the evolution for a two-qubit system. We 
summarize and discuss our conclusions in the last section. 
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II. COHERENCE MEASURES AND DYNAMIC 
EVOLUTION OF A TWO-LEVEL ATOM SYSTEM 

In this section we recall the methods to measure coher¬ 
ence in the reference basis which is due to the off-diagonal 
elements of a density matrix p , for instance, the intuitive 
1 1 norm and the relative entropy of coherence [l^ . which 
can be given by 


( 1 ) 

and 

Cre(p) = ^(Pdiag) - S{p) , (2) 

respectively. In Eqs. 0 and p indicates a state 
arbitrarily, S{p) = — Tr(/9log/9) is the von Neumann en¬ 
tropy, and Pdiag donates a density matrix that delete all 
off-diagonal elements of p. 

We will discuss the evolution of an open quantum sys¬ 
tem. The total Hamiltonian of the system with a two- 
level atom and the fluctuating vacuum electromagnetic 
field is = Hs + Hf + iJ/, where Hg and Hf denote 
the Hamiltonian of the atom and the free electromag¬ 
netic field, respectively, and Hj represents their interac¬ 
tion. Here we take a two-level atom with Hamiltonian 
Hg = ^huiQas, where ujq is the energy level spacing, and 
CTs is the Pauli matrix. The Hamiltonian representing the 
interaction between atom and electromagnetic field is 

Hi{t) =-er , (3) 

where e is the electron electric charge, er is the atomic 
electric dipole moment, E(a;) denotes the electric field 
strength. The initial state of the whole system is given 
by the total density matrix ptot = p(0) (8) |0)(0|, where 
/9(0) is the reduced density matrix of the two-level atom, 
and |0) is the vacuum state of the field. For the total 
system, its equation of motion in Schrodinger picture is 

dptotir) i 

—^— = -^[H,PtotiT)\ , ( 4 ) 

where r is the proper time. In the limit of weak coupling, 
the evolution of the reduced density matrix p(r) can be 
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written in the Kossakowski-Lindblad form 0,1111 

^ = -l[Heff,p{r)]+C[pir)], (5) 

where Hgff is the effective Hamiltonian 
1 ft t 

Heff = - /C(wo)]}cr3 , 

( 6 ) 


and 


1 X V 

C[p] = 2 X! Po-jPCTi - o-iCTjp - pcTiCTj] . (7) 

In Eqs. and ©, H is the effective energy level-spacing 
of the atom and atj is the coefficients of the Kossakowski 
matrix which can be expressed as 

O’ij — ASij , (^) 


with 

A =-[0{oJo)+ G{—i^o)] ^ B = - [0 {u}o) — G{—ujq)] , (9) 

where G{X) and /C(A) represent Fourier and Hilbert trans¬ 
forms respectively, are defined as follows: 


/ OO 

dAre*^^^G+(Ar) 

-OO 

/C(A) = — [ . (10) 

TTl J_oo W - A 

In integral the function G“''(At) is given by G~^{x — x') = 

point correlation function for electromagnetic field. Here 
|-|-), |—) denote the excited state and ground state of the 
atom. 

According to Eq. ([5]), and assuming a single two-level 
atom with initial state |V'(0)) = cos ||1)-l-e*'^ sin ||0), its 
time-dependent reduced density matrix is found to be 


, X _ 1 1 + cos^e--*^^ - 1(1 - 

PUj - 2 V ^-2Ar+^iQT+4>) 


e-2Ar-i(Qr+^) 


( 11 ) 


III. INFLUENCE OF VACUUM 
FLUCTUATIONS ON QUANTUM COHERENCE 
FOR A SINGLE-QUBIT 


li norm and relative entropy of coherence can be frozen 
during the dynamic evolution of the atom. We can do it 
by using the differential of the measures on the evolved 


Let us now analyze how the initial state and vacuum 
fluctuations affect the quantum coherence such that the 
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state with the noise parameter q that is equal to zero 


0 


and the spontaneous emission rate is given by 


5gC'(/3(<?)) = 0 , Vg e [0,1], (12) 

where < 7=1 — 6 “'’'*, t represents time and 7 is the de¬ 
coherence rate [l5l| . In this work we will consider two 
situations, where the atom interacts with a electromag¬ 
netic field without a boundary and a electromagnetic field 
with a boundary, respectively. 


7o = 


3 

E 


iTrenhc^ 


(15) 


Then we obtain the coefficients of the Kossakowski ma¬ 
trix Uij and the effective level spacing of the atom respec¬ 
tively 


A. In the case without a reflecting boundary 


^( 0 ) ^ 5 ( 0 ) ^ ^ 


In order to calculate the spontaneous emission rates of 
a static polarizable two-level atom in the electromagnetic 
field without a boundary, we consider the trajectory of 
the atom 


t(r) = T , x{t) = xo , ?/(t) = yo , z{t) = zq ■ (13) 

According to the trajectory of the atom and using the 
two-point function of the electromagnetic field in the un¬ 
bounded space [ 2 ^, we can easily obtain the field corre¬ 
lation function in the frame of the atom 


G+{x-x') = Y^ 
2=1 


e^c|(-|r.|+)p 
hn'^eoicAT — ie)'^ 


(14) 


J 


fl(0) = Wo + 


70 p 
2TrujQ 



UJ +UJo 


1 

W — Wo 


a7) 

where the superscript 0 indicates the vacuum fluctuations 
in the unbounded space. 

Using Eqs. (HH), (uni) and (HU, the 1 1 norm and relative 
entropy of coherence for the initial single-qubit are found 
to be 


Q^(p)(o) = Isinee-^^o^l , (18) 


and 


C're(p)^°^ = -i[l-bcos6»(l - g) - g]log 2 ^[1-bcos6l(l - g) - g] 

- ^[1 - cos 6»(1 - g) -b g] log 2 ^[1 - cos 6{l - q) + q] 

-b i(l-bVI)log2i(l-b VI)-bi(l-\/I)log2i(l-yi) , (19) 

where g = 1 — , and L = (1 — cos^ 0)(1 — g) -b [cos 0(1 — g) — g]^. _ 


Now let us analyze under which dynamical condi¬ 
tions the quantum coherence of the single-qubit system 
is totally unaffected by noise. By applying Eq. dUD to 
Eq. (HH), the g derivative of the lx norm can be calculated 
as follow 


I sin 0 | 


( 20 ) 


As the li norm of coherence is dependent of off-diagonal 
elements of a density matrix, we come to a conclusion 
that the necessary and sufficient freezing condition for h 
norm under the dynamic evolution of the atom is sin 9 — 
0 in the initial state. 

Similarly, substituting Eq. (1191) into Eq. © we have 


1 — cos 0(1 — g) -b g 
1 -b cos 0(1 — g) — g 


2(1 -b cos0)“1-\/L 


= 0 . 
( 21 ) 


Erom which we find that such a measure is frozen through 
the dynamic evolution only when cos 0 = 1 or cos 0 = — 1 
(trivial because the initial state is incoherent). 


Above discussions tell us that under the dynamic evo¬ 
lution leave coherence invariant only when the initial 
state is incoherent. In other words, the dynamics of 
quantum coherence can not be frozen in the studied open 
quantum system, which is due to the fact that the coher¬ 
ence is affected by the interaction between the two-level 
atom and the fluctuating vacuum. Unlike the behavior 
of a single-qubit state [l^, we find that the h norm of 
coherence can be frozen for specific initial states, while 
there is no nontrivial freezing of the relative entropy of 
coherence. 
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B. In the case with the presence of a reflecting 
boundary 

Now, we calculate the quantum coherence of the two- 
level atom in the vacuum fluctuations which existing a re¬ 
flecting boundary. To do so, we consider the case that the 
atom is placed at a distance zo from the boundary [^ . 


Then the coefficients of the Kossakowski matrix a^- and 
the effective level spacing of the atom can be calculated 
as 

A = B = ^ll-J2a^M^o,zo)], (22) 


J 


/•CX3 / I 


CdO Cd — (jJo 


(23) 


with 


fxi^jZo) — fy(X, Zq) — 


3c3 

IGA^Zq l 


2Azo 2Azo f 4A^Zn 
-cos 


c V 


° - 1 1 sin 


2Azo 


3c^ 


'0 L 


2Xzo 2Xzo . 2Xzo 

-cos-sin- 

c c c 


(24) 


r 


where, tti = |(—|ri|-|-)p/|(—|r|-|-)p represents the relative 
polarizability, which satisfies /i(a;o,^o) 

are oscillating functions of distance zq with a position- 
dependent amplitude. Hereafter, for simplify we abbre¬ 
viate <^i/i(‘j-'0j Zq) as /(wq, Zq). Then the h norm with 
the presence of a reflecting boundary can be calculated 
using Eq. (HD 

C7q(p) = |sin6»e-^'^“d-/(‘^o,^o)]r| ^ ^25) 

By applying Eq. (l25ll to Eq. (USD, we obtain the q deriva¬ 
tive of the 1 1 norm of coherence 

|sin0|[l-/(a;o,Zo)](l-9)-5[i+/(“o.-o)] = 0. (26) 

Comparing Eq. (l26t with (1201) . we can see that the term 
I —f(uJo, Zo) in Eq. (IMl) in fact decides the dynamical con¬ 
ditions so that the quantum coherence of the single-qubit 
system is totally unaffected by noise with the presence 
of a reflecting boundary. Let us now first analyze the 
limit situation when the atom is placed very far from the 
boundary, which leads to fi{uJo,zo) —?> 0, and makes this 
case reduces to that of the unbounded Minkowski vac¬ 
uum. As the polarizations of the atom are in different 
directions and behave differently. To show the proper¬ 
ties of the h norm, we plot it as the function of q and 
Zo for (ax,ay,a^) = (1, 0,0), (0,0,1), (1/3,1/3,1/3) in 
Fig. (HD, which corresponding to the parallel, vertical, 
isotropic polarization cases respectively. From Fig. (1) 
we can see that due to the presence of the boundary, the 
vacuum fluctuations is subjected to influenced, and the h 
norm of quantum coherence can be frozen under suitable 
condition. 

Case (i) It is worthy to notice that no matter what 
direction is the atom polarized in the xy plane, the re¬ 
lations az = 0 and zq —>■ 0 can be satisfied, which lead 



FIG. 1: (color online). The h norm of coher¬ 

ence as the function of q and zo for [ax,ciy,az) = 
(1,0, 0), (0, 0,1), (1/3,1/3,1/3) respectively. We set sinS = 1 
and the values of all the curves are in the unit of c/ujq. 

to fxiu}o,zo) = fy{u}o,zo) = 1 in Eq. ([SSD. Then Ci^(p) 
is independent of the environment, which signifies that 
the h norm of coherence for the single-qubit system will 
be frozen when the parallel polarized atom close to the 
boundary, protecting the important quantum resources 
of coherence effectively. 

Case (a) If the atom is polarized in z axis direction, 
we cam get ax = oty =0, which means such a measure 
is frozen only when trivially the initial state is incoher¬ 
ent, while the decoherence rate decreases slower and the 
quantum coherence get protected by the presence of the 
boundary as compared to the unbounded case. 

Case (a) If the polarization is in an isotropic polar¬ 
ization direction (ax = ay = az — only trivially 
when the initial state is already incoherent leads to the 
coherence be frozen, otherwise the quantum coherence 
become protected in some degree due to the presence of 
the boundary. 
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In the same way, the relative entropy of coherence with Eq. m 
the presence of a reflecting boundary will be found using 


Cre(p) = - 

1 

2 


1 + cos 6 (1 — q') — q' 
1 — COS 9 (1 — q') + q' 


logs 


logs 


1 + COS 6 (1 — q') — q' 
1 — COS 9 (1 — q') + q' 


+ 2 ^^ y/M) log 2 —(1 + yfM) + —(1 — y/M) log 2 —(1 — y/M) 


(27) 


where ^'(t) = 1 — e 7 = 7 o[l — f{uJo,zo)], and M = Substituting Eq. (071) into Eq. (fT7|) . we obtain the q 
(1 — cos^ 9){\ — q') + [cos 0(1 — q') — q']'^ respectively. derivative of the relative entropy of coherence 


J 


[i-/(wo,zo)](i-g)' 


- cos - 

r 


(28) 



FIG. 2: (color online). The relative entropy of co¬ 

herence is the function of q and zo for {ax,ay,az) = 
(1,0,0), (0,0,1), (1/3,1/3,1/3) respectively. We set sin0 = 1 
and the values of all the curves are in the unit of c/ujq. 

Let us consider the properties of the relative entropy 
of coherence with the presence of a reflecting boundary. 
We plot it in Fig. (ED, which is the function of q and zq 
for (a,,, Qfy, a^) = (1,0, 0), (0,0,1), (1/3,1/3,1/3), corre¬ 
sponding to the parallel, vertical, isotropic polarization 
cases, respectively. Similar to the h norm, no matter 
what direction is the atom polarized in the xy plane, the 
relative entropy of coherence for the single-qubit system 
will be frozen when the parallel polarized atom extremely 


close to the boundary. However, the relative entropy 
of coherence will be protected in other situation. In a 
word, due to the modification of the vacuum fluctuations 
caused by the reflecting boundary, the vacuum fluctua¬ 
tions induced the quantum coherence becomes position 
dependent. 

IV. INFLUENCE OF VACUUM 
FLUCTUATIONS ON QUANTUM COHERENCE 
FOR THE TWO-QUBIT SYSTEM 


Now we consider a two-qubit initial state with max¬ 
imally mixed marginals states, also known as Bell- 
diagonal {BD) states The density matrix can 

be represented as po = erf), where 

{erf} and {erf} denote the Pauli matrices. The initial 
state can also be expressed by the vector c = ( 01 , 02 , 03 ), 
which is 



/ 1 + 03 
0 
0 

V Oi - 02 


0 

1 - 03 
Ol + 02 
0 


0 

Ol + 02 
1 - 03 
0 



(29) 


According to Eq. ([SD, the corresponding time-dependent 
density matrix is 


J 


Pab{t) = -{af 0 erf -I- ciO 


-2At 


^-2At , 


COs(QT)crf ® (jf + 020 

A ^ „B I „ „-2At 


-2At 


cos{flT)af 0 erf -I- 036 ^"^^erf (g) er. 


—020 sin(QT)erf (g) af-I-OlO sin(f2r)erf 0 erf — ^(I — 0 

./I 


^fl- 


(30) 
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Comparing to the analysis of the single-qubit system 
above, we have to consider the influence of the initial 
state and vacuum fluctuations on quantum coherence of 
the two-qubit system. We also interested in under which 
dynamical conditions the coherence is totally unaffected 
by noise for the two-qubit system. In the following we 
adopt a bona fide distance-based measure of quantum 
coherence 0 

CD{p{q))=CD{po) , (31) 

which is frozen forever for any q £ [0,1], or equivalently 
frozen for any r [H, meaning that the quantum coher¬ 


ence of the system is independent of evolution. 

Substituting Eq. (1291) into Eq. ([T]) and Eq. ([2]), the h 
norm and relative entropy of the initial two-qubit state 
can be easily obtained as 


C';i(Po) = ^(|ci-F c 2|-h |ci - C2I) , (32) 


and 


Cre(Po) = - [ 2 (l-FC 3 )l 0 g 2 i(l-f C 3 )-h 2 (l-C 3 )l 0 g 2 i(l-C 3 )] 


+ [1 + C3 + (ci - C2)] log2 i[l -h C3 -h (ci - C2)] -h [1 -F C3 - (ci - C2)] log2 i[l -h C3 - (ci - C2)] 
+ [1 “ C3 -|- (ci -|- C2)] log2 -[1 — C3 -I- (ci -I- C2)] -I- [1 — C3 — (ci -I- C2)] log2 -[1 — C3 — (ci -|- C2)] 


(33) 


A. In the case without a reflecting boundary 

Let us consider the case of the vacuum fluctuations 
without a boundary. Noting that the electric-field two- 
point functions which without a boundary and the 
trajectory Eq. (USD are considered above. As a result, 
the h norm and relative entropy under the evolution are 


found to be 

Ci^{pab)^°^ = -h C 2 I -h |ci - C 2 I) . (34) 

and 


Cre{pab)^°'^ = i{ “ - 9) - 9] log2 ^[1 + C3(l - g) - g] - [1 -H 03(1 - q) + q] log2 ^[1 -t 03(1 - q) + q] 

- [1 - C3(l - g) - g] log2 i[l - C3(l - g) - g] - [1 - 03(1 - g) -h g] log2 ^[1 - 03(1 - q) + q] 

-h (1 -h C3(l -q)- N) log2 i [1 -h 03(1 -q)-N]+{l + 03(1 - q) + N) log2 ^ [l -F 03(1 - q) + N] 

-h (1 - C 3 (l -q)- N) log2 i [1 - C 3 (l -q)-N]+{l- 03(1 - q) + N) log2 ^ [l - 03(1 - g) -f A] | 


(35) 


where N = y^(l — g)(ci -I- 02)^ -I- g^, and the superscript 
0 indicates the vacuum fluctuations in the unbounded 
space. Now let us analyze under which dynamical condi¬ 
tions the quantum coherence of the two-qubit system is 
totally unaffected by noise. By applying Eqs. (l32ll — (l35ll 
to Eq. we find that the coherence oi BD states will 
be frozen for any bona fide distance-based measure only 


when 

Cl = C2 = 0 . (36) 

Therefore, there is no exist trivial freezing of coherence 
for the dynamics of the two-qubit system, which means 
that the quantum coherence of BD states must be af¬ 
fected by noise with time. 
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B. In the case with the presence of a reflecting 
boundary 

However, since vacuum fluctuations will be modified 
for the present of boundary in the vacuum, it is neces¬ 
sary to examine how the vacuum fluctuations affect the 
quantum coherence. According to the electric-field two- 
point functions with the presence of a boundary [2^ and 

I 


the trajectory Eq. m, the 1 1 norm and relative entropy 
of the coherence under the evolution can be easily ob¬ 
tained 

CiApab) = ie-^[i-/(“«’^«)l"(|ci+C2| + |ci-C2|) , (37) 
and 


Cre{pab) = + C3(l - q') - q'] log 2 ^[1 + C3(l - q') - q'] - [l + C3(l - q') + q'] loga ^[1 + C3(l - q') + q'] 

- [l - C3(l - q') - q'] log2 i[l - C3(l - q') - q'] - [l - €3(1 - q') + q'] loga ^[1 - 03(1 - q') + q'] 

+ (1 + C3(l — q') — S) log2 ^ [1 + ^3(1 — q') — S'] -I- (1 -I- 03(1 — q') + S) log2 ^ [l + *^3(1 — q') + S] 

-f (1 - 03(1 - q') - S) log2 i [1 - €3(1 - q') - S] -f (1 - C3(l - q') + S) log2 | [l - 03(1 - 9') + -S'] | , 

(38) 


where S = y^(l — q'){ci + cAA + q'^- According to 
Eq. dSIl), the dynamical condition of the quantum co¬ 
herence of the two-qubit system is determined by the 
term 1 — f{uJo,zo) when a boundary is present. Simi¬ 
lar to the analysis of the case of the single-qubit system, 
we can see that whatever direction the polarization is 
in the xy plane one has = 0, resulting in the factor 
1 — /(wo,zo) = 0. This dynamical conditions for freez¬ 
ing the quantum coherence can be fulfilled if the atom 
is transversely polarized and very close to the boundary, 
which means that the coherence of a two-qubit system 
can be frozen. Otherwise, the coherence must be affected 
by noise under the dynamic evolution in open quantum 
systems and can only be protected in some degree. 

V. CONCLUSION 

In the framework of open quantum systems, we have 
investigated the dynamic evolution of a single-qubit sys¬ 
tem and a two-qubit system when the atom is cou¬ 
pling with the bath of fluctuating vacuum electromag¬ 
netic field. Two different situations, the atom interacts 
with an electromagnetic field without or with the pres¬ 
ence of a reflecting boundary, have been studied. For a 
single-qubit system, we find that for the case of without a 
boundary, the quantum coherence will be destroyed due 


to the interaction between the atom and the electromag¬ 
netic field. However, with the presence of a boundary, 
the quantum coherence becomes dependent on position 
and atomic polarization. The dynamical conditions for 
the insusceptible of coherence are fulfilled only when the 
atom is close to the boundary and is transversely po¬ 
larizable. Otherwise, the quantum coherence can only 
be protected in some degree in other polarizable direc¬ 
tion. Similarly, for a two-qubit system, the vacuum fluc¬ 
tuations always affect the quantum coherence without a 
boundary. With the presence of a boundary, the quan¬ 
tum coherence of two-qubit states with maximally mixed 
marginals is shielded from the influence of the vacuum 
fluctuations of the electromagnetic field when the atom 
is close to the boundary and is transversely polarizable. 
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